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^ I Abstract. Freeze-out of particles in relativistic hydrodynamics is considered across 

^ ■ a 3- dimensional space-time hypersurface. The conservation laws for time-like parts of the 

O . freeze-out hypersurface require different values of temperature, baryonic chemical potential 

Q>^ , and flow velocity in the fluid and in the final particle spectra. We analyze this freeze-out 

discontinuity and its connection to the shock-wave phenomena in relativistic hydrodynamics. 
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1. Introduction 

Relativistic hydrodynamical model |]I|] has been widely discussed in recent years within their 
connection to high energy nucleus- nucleus (A+A) collisions (see, for example, 0). This 
approach is the most straight and unique way to investigate the hadron-matter equation 
of state (EoS) and possible phase transition into a deconfined state called the quark-gluon 
plasma. The aim of the relativistic hydrodynamical model of A+A collisions is to elucidate 
the properties of the highly excited matter comparing the model results with experimental 
observables. The system evolution in relativistic hydrodynamics is governed by the energy- 
momentum tensor 

T^- = {e + p)uf'u'' - pgf"" (1) 

and conserved charge currents. The baryonic current, 

f = nu^ , (2) 

plays the main role in the applications to A+A collisions. The hydrodynamical description 
includes the local thermodynamical variables (energy density e, pressure p, baryonic density 
n) and the collective four-velocity m'^ = (1 — v^)~^/^(l, v). The continuous flows are the 
solutions of the hydrodynamical equations 

d.T^'^ = , d,f = , (3) 

with specified initial conditions. To complete the system of hydrodynamical equations (|^) 
the fluid EoS has to be added as an input p = p{e,n). It is usually used in a more con- 
venient form oi p = p{T,^b) with temperature T and baryonic chemical potential /i^ as 
independent variables. Other thermodynamical functions can be found from the following 
thermodynamical identities (s is the entropy density): 

n = - — , e = Ts + fiBn -p ■ (4) 

\ohbJ T 

Besides the fluid EoS and the initial conditions for the fluid evolution the third crucial 
ingredient of the hydrodynamical model of A+A collisions is the so-called freeze-out (FO) 
procedure, i.e. the prescription for the calculations of the final hadron observables. Parti- 
cles which leave the system and reach the detectors are considered via FO scheme, where 
the frozen-out particles are formed on a 3-dimensional hj^ersurface in space-time. The 
FO description has a straight influence on particle momentum spectra and therefore on all 
hadron observables. During last three years an essential progress has been achieved in the 
formulation of the FO procedure in relativistic hydrodynamics. A generalization of the well 
known Cooper-Frye (CF) formula P| to the case of time-like FO hypersurface was suggested 
in Ref. 0]. The new formula does not contain negative particle number contributions on 
time-like FO hypersurface appeared in the CF procedure from those particles which cannot 
leave the fluid during its hydrodynamical expansion. The FO procedure of Ref. [Q has been 
further developed in a series of publications [5-10]. Ref. [8] presents also the kinetic approach 
to the FO problem. 

There is a new distinct and unusual feature of the proposed FO procedure. The param- 
eters of the distribution function for the final particles should be different from those in the 




fluid to satisfy tlie conservation laws. The particle emission from the time-like parts of the 
FO hypersurface looks as a 'discontinuity' in the hydrodynamic motion. We call it the FO 
shock. In the present paper we analyze a connection of this new type of discontinuities to a 
general problem of shock waves in relativistic hydrodynamics. 

2. Freeze-out procedure 

In the hydrodynamical model of A+A collisions the fluid expansion has to be ended by 
the correct transformation of the fluid into free streaming hadrons. The most important 
requirement for this FO procedure is to satisfy energy-momentum and charge conservations 
during the fluid transition into free final particles. It is assumed that there is a narrow space- 
time region where the mean free path of the fluid constituents increases rapidly and becomes 
comparable with the characteristic size of the system. The local thermal equilibrium is 
supposedly maintained in this intermediate region. In practice, one considers a "zero width" 
approximation and introduces a FO surface, so that particle distributions remain frozen-out 
from there on. 

In his original paper Landau []I| defined the FO hypersurface by the condition T(t, x) = T*, 
where t and x are the time and the space coordinate respectively, and T* is the fluid FO 
temperature chosen to be approximately equal to the pion mass. The FO procedure, first 
introduced by Milekhin |]Tl| , was based on the intuitive physical assumption of isotropic ther- 



mal particle spectra in the rest frames of each fluid element. The final hadron observables 
was then obtained by the Lorentz transformations of these spectra with known hydrody- 
namical velocities and integration over all fluid elements. The FO procedure was improved 
further by Cooper and Frye [Q and this method has been used ever since. Final spectrum 
for z-th type hadron is expressed by formula 0: 

where do"^ is the external normal 4- vector to the FO hypersurface Sj, 0j denotes the local 
thermal Bose or Fermi distributions, di is the degeneracy factor for particle i with the chem- 
ical potential /Xj. The particle spectrum (^ is known as the Cooper-Frye (CF) distribution 
function. The chemical potential fii is equal to 

/Wi = bifiB + Sifis + QifJ'Q , (6) 

where 5j, Sj, qi are the baryonic number, strangeness and electric charge of i-th particle. 
The strange chemical potential fis and electric chemical potential yUg are always chosen to 
make the total strangeness to be equal to zero and the ratio of electric charge to the baryonic 
number to be fixed at their initial value given by A+A collision. Therefore, only the baryonic 
chemical potential fiB is considered as independent thermodynamical variable. 

The CF formula (|^) corrected the Milekhin's approach. The matter is that the external 
normal vector da^ equal to {—dx, dt) in 1+1 dimensions does not in general become equal 
to da* = {—dx*,0) after the Lorentz transformation into the fluid element rest frame. The 
non- simultaneous particle FO in the rest frame of the fluid element modifies the thermal 
distribution function and causes a nonzero momentum of the fluid element in its rest frame 
("local anisotropy", see Ref. |T2[ for details). 



The initial conditions of the hydrodynamical motion are given on the initial hypersurface 
Sj„. The final hypersurface Sj should be closed to Sj„ and in general Sj consists of both 
space-like (s.l.) and time-like (t.l.) parts. Using the Gauss theorem and Eq. (|]) one can 
easily prove the energy-momentum and charge number conservations with CF FO formula: 



(27r)3 hf "^ J ^^y T* J ' 

where bi is the baryonic number of i-th particle and the summation over all particles and 
the corresponding antiparticles is assumed in the right hand side of Eqs. (0,|). Note that 
fiuid EoS at the FO is assumed to be an ideal gas EoS for a mixture of different particle 
species. The hydrodynamical quantities T^'^ (|1]) and j^ (^ are identically presented on the 
FO surface Sj in terms of final particle distributions (pi according to Eqs. (|^,|). 

The CF formula (^ still does not provide a complete solution of the FO problem. The 
FO surface consists of t.l. parts too. Eq. (§) can not however be used for a t.l. FO surface 
(i.e., s.l. normal vector da,): free final particles "return" to the fiuid if da,k^ < 0, and this 
causes unphysical negative contributions to the number of final particles. A first attempt to 



improve the CF formula for t.l. FO hypersurface was done by Sinyukov in Ref. |I3|. The 
modified FO procedure and new formula for the final particle spectra emitted from the t.l. 
FO hypersurface was proposed in Ref. @]: 

k^"-^ = A / da,k^ ^^ f ^^^^1 e(da^kn . (9) 



dk^ (27r)3 hj M ^* \^ r^^ 



a 



Eq. (^ looks like CF formula (|^), but without negative particle numbers that appear in (|^) 
for t.l. FO surfaces. These negative contributions are cut-off by the ^-function in Eq. @. 
We'll call Eq. (0) the cut-off (CO) distribution in what follows. An inclusion of the CO FO 



into the self-consistent hydrodynamical scheme is considered in Ref. [|10|- The distribution 
function for the free particles in Eq. (|^) contains the new parameters Tg, nf, Ug. We'll briefly 
call the free particle state as the 'gas' in order to distinguish it from the ordinary fluid. 
Note, however, that flnal particle system differs from the normal fluid and gas. It has the 
non-thermal CO distribution function which is frozen-out as all further particle rescatterings 
are 'forbidden' in the post-FO gas state. 

The conservation laws across FO hypersurface S/ can be written now in their local form 
as 

(27r)3 ^^ J k- 

E^^-. /^^"^0^(^y^| (11) 




E jA-.dcr, I ^ k^k^ 0. ( I^r^^ ) e{da^k'^) 



These conservation laws were derived together with the equation of motion of the fluid in 
Ref. []10[ and should be solved simultaneously with the latter ones in order to find the FO 



hypersurface in the practical hydrodynamical calculations. 

The presence of the 6'-function in the right hand side of Eqs. ( p^miTlD leads to the dis- 



continuity between the fluid T/, /i/, Mj and the gas Tg, ^f, u'^ variables to satisfy the energy- 
momentum and baryonic number conservation. We call this discontinuity the FO shock. 
Note that the strangeness flow across the FO surface is equal to zero and electric charge flow 
is proportional to the baryonic number flow. Therefore, the conservation of strangeness and 
electric charge is also fulfilled. 

Another FO discontinuity problem occurs in the hydrodynamical model calculations when 
the EoS is different on the two sides of the FO front (for the final gas of free particles one 
must use the ideal gas distributions). The conservation laws on the FO surface in this case 
should lead to the discontinuities between fluid and gas of final particles across both t.l. and 



s.l. parts of Tif ^M. The discontinuity across the s.l. part of Sj corresponds to the so-called 



time-like shocks |JT^ (i.e., (icr^ is t.l. vector) and will be also discussed in the next Sections. 

3. Relativistic shock waves 

Equations (^ are the differential form of the energy-momentum and baryonic number con- 
servation laws. Along with these continuous flows, the conservation laws can also be realized 
in the form of discontinuous hydrodynamical flows which are called shock waves and satisfy 
the following equations: 

Tr A, = T^'A, , no<A^ = nu^A^ , (12) 

where A^ is the unit 4-vector normal to the discontinuity hypersurface. In Eq. (|12|) the zero 
subscript corresponds to the initial state ahead of the shock front and the quantities without 
an index are the final state values behind it. 

The important constraint on the transitions (|12]) is the requirement of non- decreasing 
entropy (thermodynamical stability condition): 

su^A^ > SoKA^ . (13) 

We consider one-dimensional hydrodynamical motion in what follows. In its usual sense 
the theory of the shock waves corresponds to the discontinuities on the surface with space- 
like A^, i.e., the shock-front velocity is smaller than 1. In this case one can always choose the 
Lorentz frame where the shock front is at rest. Then the hypersurface of shock discontinuity 
is Xsh = const and A^ = (0, 1). The shock equations (|12]) become 

T^' = T'\ Tl' = T" , ny^ = nu' . (14) 

From Eq. (p!^ one obtains 

^2 ^ (p-Po)(g+Po) ^2 ^ (p-Po)(£o+p) ,^^. 



Substituting (|T5]) into the last equation in (0) we obtain the well known Taub adiabate 
(TA) equation p[ 

n'X^ - nlXl - (p - Po)(X + Xo) = , (16) 

in which X = (e + p)/n'^, and TA therefore contains only the therniodynamical variables. 
Eq. ([T6|) can be written also in the form 



2 _ 2 {^ + Po){.e + p) 



{so+Po){eo+P) 



If the equation of state p = p{e,n) is given, then Eq. (|T^ defines the function n = riT^s) 
which is the TA (|T6|) in the {e,n) plane. The TA in the {s,p) plane is defined by the 
function Pt{^) = pi^i^Ti^))- In the case of zero baryonic number it is just reduced to the 
EoS Pt{^) = p{^)- The point (^ojPoj'^o) is called the center of the TA. The shock transition 
from the state {eo,Po,no) to {ei,pi,ni) is mechanically stable if for all e between ^o and ei 



the following inequality is valid |T7|: 



(ei-Eo) [Pcr{e) - Prie)] > 0, 
where 



with 



Pcrie) = A - -^ (19) 

A = ''^' 7 ^°^° , C={A + e,){A-p,) = {A + e,){A-p^). (20) 

ei-Eo+Po- Pi 

Therniodynamical stability (0) follows from the mechanical stability (|15), and the inverse 
statement is not in general true |l^. It follows from (|18D that the necessary and sufficient 



condition for the mechanical stability of shocks is that the compression {ei > Eq) TA pris) 
goes under the critical curve Pcri^) ([T9|) and the rarefaction {ei < Eq) TA goes above Pcri^) 
for all e between the initial ^o and final ei points. The consequences of these relations are the 
well-known inequalities for the speed of sound, Cg = [{dp/de)s/nY^'^, and the fiow velocities 
at both sides of the shock front in its rest frame 

CsO < Vo , Csl > Vi , (21) 

which are the necessary conditions of the general mechanical stability criterion (|TB|). 
The function 

defines the thermodynamic properties, i.e., the medium is said to be thermodynamically 
normal if the quantity zt is positive, and thermodynamically anomalous when zt < 0. In 
the normal case zt > the compression shocks are stable whereas the rarefaction shocks 
become stable in the thermodynamically anomalous media with zt < 0. Note that the 
function Pcri^) (|T9|) is the solution of the equation ZT^e) = 0. The case of varying sign of 



zt{s) which is typical for the phase transitions are discussed in Refs. |T9|, ^ . 

Let us consider now the discontinuities on a hypersurface with a t.l. normal vector A^. 
This possibility was suggested by Csernai in Ref. |T^. We call them t.l. shocks. In this case 



one can always choose another convenient Lorentz frame ('simultaneous system') where the 



hypersurface of the discontinuity is tgh = const and A^ = (1, 0). Equations ([T2|) become then 

T°° = T°°, T„i° = T^o , ny, = nu"^ . (23) 

From Eq. (p3D we find 



„2 {e-eo){eo + p) ^2 (e - eo)ie + Po) 



° {p-Po){e + po) ' (p-Po)(eo + p) 



where we use " ~ " sign to distinguish the t.l. shock case (p^) from the standard s.l. shocks 
(p!5|). Substituting (|2^) into the last equation in (^) one finds the equation for t.l. shocks 
which is identical to the TA of Eq. ([T6|). We stress, however, that the intermediate steps are 
quite different. The two solutions, Eqs. ( ^41) and ([15|) , are connected to each other by simple 
relations, 

vl = \, ^' = 4 ' (25) 

yZ yZ 

between the velocities of s.l. and t.l. shocks. These relations show that only one type of 
transitions can be realized for the given initial and final states. Physical regions [0, 1) for 



■^oi 



([15|) and for v'^,v'^ (|2^) can be easily found in the (£:-p)-plane shown in Fig. 1. 



If one takes as the initial and final states only states which are thermodynamically equi- 
librated, the TA passes then through the center of TA {eo,po)- For any physical EoS with 
< Cs < 1 TA lies as a whole in the regions I and IV in Fig. 1 |]18[, i.e. only compression 
and/or rarefaction s.l. shocks (with s.l. normal vector A^) are permitable. 

The only way to make the t.l. shocks to be possible is to allow the metastable initial 
and/or final states. The TA then no longer passes through the initial point (£o,Po) and new 
possibilities of t.l. shock transitions (^,0) to regions III and VI in Fig. 1 would be realized 
(see, e.g., the t.l. shock hadronization of supecooled quark-gluon plasma in Ref. pT[|). 

4. Freeze out shocks 

The conservation laws ( [TI1|JTT| ) between the fluid and the free streaming particles are sim- 
plified across the s.l. FO hypersurface. In this case the ^-function in the right hand side 
of Eqs. (p!0| , pTD is equal to unity for all momenta of final particles. The CO formula (P) is 



reduced to the CF distribution (^) at s.l. parts of the FO hypersurface. It was shown in the 
previous Section that t.l. shocks (shocks across s.l. hypersurface) are not permitable for any 
physical EoS, in particular for the ideal gas EoS for the fluid and the gas of free particles. 
It means that the only solution of Eqs. (|lO| , pTD is 



Tf = Tg, ^4 = ^4, u^f = u^g. (26) 

No discontinuities between the fluid and the gas parameters can exist, and the particle 
emission from s.l. hypersurface is a smooth FO according to the CF formula (|^). We do not 
discuss here the possibilities when the fluid FO takes place from the state which is out of 
local thermodynamical equilibrium. 

For the t.l. parts of the FO hypersurface Sj a rather different solution should be realized. 
The ^-function in Eqs. (|TU|,|TT]) equals to zero when da^k'^ is negative. Therefore the smooth 



solution ( p6D for Eqs. ( [TO| . [Tl| ) does not exist. In this case the only possibility to satisfy the 
conservation laws (^^,0) is to introduce a discontinuity between the fluid Tj, fi{ , u'^ and gas 
Tg, fif, u'i parameters. Let us consider this FO shock on t.l. hypersurface in more detail. To 
obtain the analytical solution of the problem we restrict ourself to the case of zero baryonic 
number n = and consider massless particles with the Jiitner distribution function 

, (^) ^ exp (- ^) . (27) 

Similar problem, but for the conserved particle number (not the charge!), was considered in 
Ref. i. 

Substituting it into both sides of the Eq. (|lOl), one finds that the energy-momentum 
conservation between the fluid and free particles along the t.l. hypersurface acquires the 
form: 



da.Tg^- =da,J — k^k^ <P {-^\ e{daxk^) . (28) 



The energy-momentum tensor T'^'^ in Eq. (pHj) has the form of Eq. (|T]) with the functions 

e{Tf) = Sf and p{Tf) = pf given by 



CXD 



27r2 7o '^ " ^"^ \ f] ~ ^ 



e{T) = 3p{T) = -^ / k^dk k exp -- = — T^ . (29) 



In the rest frame of the FO front, i.e. in the Lorentz frame where (ia^ becomes equal to 
(0, dt), the Tg^ and T^^ components of the gas energy- momentum tensor in Eq. (|28D can be 
rewritten then in the form 

Tf = isl+pl)uyg, Tl' = {el+pl)u]u] +pl. (30) 

Note that such a representation can be given only for the T^-^ components which enter the 
energy-momentum conservation (PH]). This form coincides with that of Eq. ([^), but with 
effective values of £* and p* which are found to be equal to 



^;-w)^. p;^Pm '^+"'f-"»' . (31) 



where Vg is the velocity parameter of the gas in the FO shock rest frame. The functions e 
and p in the right hand side of Eqs. (|3l[) are given by Eq. (|^). 

Due to the same formal structure of T° ^, T^ ^ ( |5D| ) and T9^, TI^ (Q) one obtains the solution 
of Eq. (H) 

^ {ef-el){ef+pl)' ' (e/ -£;)(£; +p/) ' 

which is similar to Eq. (|T3|). Note, however, that in contrast to Eq. (|T3|) the values of e* and 
p*g in Eq. (|32D are not just the thermodynamical quantities, but depend also on Vg. 



It should be emphasized that Eq. (^) can also be introduced for the distribution functions 
of massive and charged particles [0 . In the latter case the energy-momentum conservation 
leads to the familiar expressions ( ^2l) for the effective energy density and pressure, and the 
charge conservation leads to the TA equation ([T6|) or ([1^) for the effective charge density. 

We fix the FO hypersurface Sj by the condition Tg = const. The analytical solution of 
Eq. (p2D can be presented then in the form 






It is instructive to compare these results with shock-wave solution (^) for the same ideal gas 
EoS (^91). We fix the temperature Tg of the final state in s.l. shock wave (0) and consider 
Vg and Tf dependence on Vf-. 




3(1 -t; 



Figs. 2 and 3 show the dependences of Vg and Tf on Vf for the fixed value of Tg in the 
final state, both for the FO shock ( ^3|) and for the normal shock wave (0). The kinematical 
restrictions on the gas velocity give the same value of the minimal fluid velocity in both 
shock transitions, (f/)mm. = Pf/^f- It equals to 1/3 for the considered ideal gas EoS (|29| ) 
for the fluid. Figs. 2 and 3 indicate that at low values of the fluid velocity, Vf < 1/^/S, 
the behavior of Vf and Tf/Tg for the FO shock (|33D and for the normal shock wave (|3^) is 
quite similar. The value of Vf = l/-\/3 corresponds to the speed of sound, c^ = 1/^/S, in the 
system with ideal gas EoS (pOl) . According to the requirements given by Eq. (|2l|) the shock 
transitions at Vf < Cg are mechanically unstable. Mechanically stable solutions at Vf > Cs 
for normal shock waves and FO shocks have qualitatively different behavior. For the stable 
normal shock wave one has Vg < Vf and Tf < Tg (see Figs. 2, 3), i.e. only compression normal 
shock wave transitions '/ '^g ' would be stable. For the FO shock transitions we find a 
completely different behavior, Vg > Vf and Tf > Tg, illustrated in Figs. 2, 3. In contrast to 
the normal shock waves, only the rarefaction FO shock transitions '/'— '■'g' ' are stable. This 
result is in agreement with an intuitive physical picture of the FO as a rarefaction process. 

The decompression degree in the FO shock can be estimated by the energy density ratio 



on both side of the discontinuity. In Appendix B of Ref. |]T0| the energy density of the gas 
was found according to the Landau-Lifshitz definition of the collective velocity [^]. In order 
to find the desired quantity it is necessary to diagonalize the energy-momentum tensor of 
the gas by the Lorentz transformation. Zero-zero component of the obtained tensor is the 
energy density. It is found to be smaller than the effective energy density e* and is given by 



T^'iv^Tg] 



/ A\2 V4 + 2A + A2 + 1 + A 



where A = 1—Vg and Vg is defined by Eq. (0). The magnitude of decompression in the FO 
shock is presented in Fig. 4. It reaches the maximal value of about 2 at the edge of of the 
stability region, Vf = l/VS. Therefore, the quantitative change of the hydrodynamical values 
in the FO shock is not negligible and it might be manifested in the various applications. 

8 



The entropy flux of tlie fluid is given by 



s'ida,, 



SfU^da 



fi ; 



(36) 



wliere Sf = [ej +pf)/Tf is tlie fluid entropy density. Similar expression is valid for the final 
'gas' state in the case of normal shock waves. The entropy production in the normal shock 
waves is given by the formula 






33/4 



^/ 



9^/ - 1 



1/4 



(37) 



shown in Fig. 5. The entropy flux of the gas with the cut-off distribution function is given 
byil: 






~t7 



In 



kuUg 



e{d(xxk^ 



(3J 



The entropy production in the FO shocks is calculated then as 



s'^da^ 
s^fda^ 



1 

2^ 



(3.;^ + 1)^(1 - 3t;;) 
(1 + Vf) 



1/4 



(39) 



The maximal entropy production in the FO shock (^) corresponds to Vj = Cg = l/"\/3 and 
can be considered as an analog of the Chapman- Jouguet point. One has to note that the 
Chapman- Jouguet point of the FO shock is also the boundary of the mechanical stability. 
It is similar to the case of the usual shocks |T^, |TB|, |T^, ^ P3[ . 

The initial {ef,pf) values in the FO shock belong to the straight line p 



plane. The respective final effective values (£q,p! 



e/3 in the e-p 
are defined by the energy-momentum con- 



servation (pOl). For the fixed FO temperature {Tg = const) the possible values of {ef,pf) and 
the corresponding values of {e*p*) are shown in Fig. 6. The point O in Fig. 6 corresponds 
to f/ — i> 1. It follows that Vg ^ 1 a.s well and the non-cut distribution function of the gas 
becomes identical to the fluid one. If the fluid velocity decreases, the gas velocity decreases 
also. The difference between the fluid and the gas becomes then more pronounced. This 
leads to a small growth of the gas effective energy (on the arc OB ) and and an essential 
growth of the fluid energy density (on the interval OA). The point A is the Chapman- Jouguet 
point, which corresponds to the maximal value of the effective energy density of the gas. The 
unstable final states also belong to the same arc BO, but they correspond to the opposite 
direction of changes. The point A is the mechanically stable initial state of the fluid with 
maximal energy density. The mechanically unstable initial states are located to the right 
hand side of the point A. 

Curve AOB in Fig. 6 is not a usual TA. The principal difference is that the gas temper- 
ature is fixed by the FO condition and there are no states with intermediate temperatures 
between the given initial and final ones. However, the mechanical stability condition ([T8| ) 
corresponds exactly to the case of usual shocks and TA. The critical curve (|T^) for points 
A and B in Fig. 6 is a tangent to the line OA at point A. The stable FO shock transitions 
satisfy the mechanical stability criterion. For the unstable FO shocks the critical curve has 



the intersection point with the hne p = e/S within the interval OA. Suppose there is such an 
intersection point D for the FO shock transition from the initial state C to the final state E of 
the arc BO in Fig. 6. Then Fig. 5 demonstrates that entropy production on the DE interval 
is smaller than that on the interval AB. Moreover, the rarefaction shock transition C— >D 
is unstable because it connects the two states of the thermodynamically normal matter. 
Therefore in the C— >-D transition the entropy decreases. Combining both results together, 
one finds that the tangent point A of the critical curve corresponds to the maximum of the 
entropy production, i.e., it is a Chapman- Jouguet point. 

5. Summary 

In the present paper the particle freeze-out in relativistic hydrodynamics was considered. For 
time-like parts of the FO hypersurface it is described by the cut-off distribution function and 
particle emission looks as a discontinuity in the hydro dynamical motion (FO shocks). The 
connection of these FO shocks to the normal shock waves in the relativistic hydrodynamics 
has been analyzed. 

The FO shock transition of the fluid to the gas of free particles with the cut-off distribu- 
tion function is mechanically stable if the fluid velocity Vf in the rest frame of the FO shock 
front is larger then the speed of sound c^ inside the fluid. This condition coincides with the 
mechanical stability criteria for the normal shock waves. In this case the particle emission 
does not influence the fluid evolution and can be selfconsitently included in the relativistic 
hydrodynamical models. We have used only the consequence (21) of the mechanical stability 
condition. Its complete formulation for the t.l. F.O. requires further studies. 

Analytical solution for the FO shocks are obtained in the case of the ideal gas EoS (P^. 
In contrast to the normal shock waves, FO shocks are found to be the rarefaction shocks: 
Tg < Tf and also p* < pf, e* < ej. For the normal shock waves such behavior corresponds 
to the matter with anomalous thermodynamical properties. The entropy increase (thermo- 
dynamical stability) takes place for the rarefaction FO shocks, but the entropy production 
(p^) appears to be quite small. 

The FO procedure across a time-like hypersurface considered in this paper could have 
rather wide field of applications. It is necessary for a correct transformation of hydrodynam- 
ical quantities into the spectra of produced particles. One example is the model of A+A 
collisions which combines hydrodynamics for the early deconfined stage of the reaction with 



a microscopic non-equilibrium model for the later hadronic stage (see Ref. p^). 
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Fig. 1. Possible final states in the e — p plane for shock transitions from the initial state 
(eo,Po)- I and IV are the physical regions for s.l. shocks, III and VI for t.l. shocks. II and 
V are unphysical regions for both types of shocks. Note, that only states with p < e are 
possible for any physical equation of state in relativistic theory. 
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Fig. 2. Gas velocity Vg as a function of the fluid velocity vj in the rest frame of the 
shock front. The dashed lines represent mechanically unstable transitions, whereas the solid 
lines show the mechanically stable FO shocks and normal shock waves. Thin lines corre- 
spond to a normal shock wave (Eq. (p^ ) and the thick ones correspond to the freeze-out 
shock (Eq. (|5BD). The circle has the coordinates (1/v^; 3 — 4/-\/3) and the square has the 
coordinates (l/-\/3; l/-\/3). 
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Fig. 3. Ratio of the temperatures on the both sides of the shock front, Tf/Tg, as a func- 
tion of the fluid velocity Vf in the rest frame of the shock front. The legend corresponds 



to the Fig. 2. The circle has the coordinates 1/a/3; (8/(3\/3))^/'^ 
coordinates (1/a/3; 1). 



and the square has the 
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Fig. 4. Ratio of the energy densities on the both sides of the shock front, as a function 
of the fluid velocity Vf in the rest frame of the freeze-out shock front. The dashed lines 
represent mechanically unstable transitions, whereas the solid lines show the mechanically 
stable FO shocks and normal shock waves. Thin lines correspond to a normal shock wave 
and the thick ones correspond to the freeze-out shock . The energy density Eg is given by 
Eq. ( |29| ) for the normal shock wave and for the gas of free particles after the FO shock it is 

defined by Eq. (p5|). 
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Fig. 5. Ratio of the entropies on the both sides of the shock front, as a function of the 
fluid velocity Vf in the rest frame of the shock front. The legend corresponds to the Fig. 
2. The upper and lower panels differ by the vertical scale. The circle corresponds to the 
maximal entropy (analog of the Chapman- Jouguet point) in the freeze-out shock and it has 



the coordinates 



1/V3; 



4(V3+1) 



and the square has the coordinates (l/\/3; 1). 
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Fig. 6. Possible initial and final states in the e—p plane for the freeze-out shock transitions. 
Initial fluid states are on the straight line p = e/3. Each initial state on OA interval 
corresponds to the one point on curve OB of the final gas states in the FO shock (see text 
also). The initial fluid states of the unstable FO shock transitions are shown by the long 
dashed line. The critical curve Pcri^) for A^B FO shock (short dashed line) is a tangent 
to the line p = e/3 at the point A which is an analog of the Chapman- Jouguet point. The 
critical curve Pcri^) for C^E FO shock (dashed-dotted line) crosses the line p = e/3 at the 
point D. 
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